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We present a systematic study of the correlation-induced corrections to the electronic band struc- 
ture of zinc-blende BN. Our investigation employs an ab initio wave-function-based local Hamil- 
tonian formalism which offers a rigorous approach to the calculation of the polarization and local 
charge redistribution effects around an extra electron or hole placed into the conduction or valence 
bands of semiconducting and insulating materials. Moreover, electron correlations beyond relax- 
ation and polarization can be readily incorporated. The electron correlation treatment is performed 
on finite clusters. In conducting our study, we make use of localized Wannier functions and em- 
bedding potentials derived explicitly from prior periodic Ifartree-Fock calculations. The on-site and 
nearest-neighbor charge relaxation bring corrections of several eV to the Hartree-Fock band gap. 
Additional corrections are caused by long-range polarization effects. In contrast, the dispersion 
of the Hartree-Fock bands is marginally affected by electron correlations. Our final result for the 
fundamental gap of zinc-blende BN compares well with that derived from soft x-ray experiments at 
the B and N K -edges. 



I. INTRODUCTION 

Electronic structure calculations for solids are com- 
monly based on solving self-consistently a single-particle 
Schrodinger equation. The solutions have the form of 
Bloch waves, i.e., they are translation symmetry adapted. 
First principles calculations of this form are success- 
fully perfo rme d at present within density functional the- 
ory (DFT^^or Hartree-Fock (HF)'^ approaches supple- 
mented by post-HF calculations. A particular object of 
interest in these studies is the correct treatment of elec- 
tron interactions in crystals with relatively weak as well 
as strong electron correlations. 

Within the framework of DFT electron correla- 
tion effects are incorporated in practice through var- 
ious exchange-correlation functionals, whereas diverse 
post-HF correlation algorithms have been designed 
to addre ss this issue with the w ave-function-based 
approach.Bllll i l 8 | 9 | io | ii | i2 | i3 | i4 | i5 | i6 | i7 | is strictly 

speaking a ground-state theory. Nevertheless the eigen- 
values of the Kohn-Sham (KS) equations have been inter- 
preted as the electronic energy bands. Recent theoreti- 
cal and computational studies on molecular systems have 
shown that in analogy with Koopmans' theorem in HF 
theory, the absolute values of the energies of the KS (spin- 
) orbitals can be indeed interpreted as approximate ver- 
tical ionization potentials. P^^ For the highest o ccupied 
KS orbital, this relation becomes exact .ESMlMIi] Yet, the 
calculated ionization potentials depend on the approxi- 
mation to the exact KS potential. 

The description of the electronic properties of many 
elemental metals and various covalent and ionic com- 
pounds by means of DFT has appeared to be very suc- 
cessful. Recent extensions of the theory to the time- 



depe ndent d omain, time-dependent DFT and current 
Qprp|22i23i24| jjg^yg g^iso allowed for a more rigorous treat- 
ment of optical excitations. By circumventing the calcu- 
lation of the many-body wave function, ground or excited 
state, DFT offers a relatively simple approach to the elec- 
tronic structure of solids. Such studies are commonly car- 
ried out within the local-density approximation (LDA^P^ 
or generalized gradient approximation (GGA).^® 

It is well known, however, that the LDA-DFT and 
GGA-DFT based calculations underestimate the band 
gaps of insulators and semiconductors. It is also known, 
that this failure is only partially related to the use of 
approximate exchange-correlation functionals. What a 
more complete determination of energy gaps does require 
is an exchange-correlation potential which is orbital de- 
pendent. The polarization and relaxation cloud around 
an extra electron placed into the conduction bands or 
a hole created in the valence bands is of different na- 
ture than that of th e elec trons in the ground state of 
the solid; see, e. g.,^^^ The addition or removal of 
an electron to or from the infinite system leaves the 
electron density unchanged and therefore, the exchange- 
correlation potential in LDA-DFT, which depends on the 
local densities only, can not d escrib e properly the corre- 
lation hole of the extra charge .'^^'^Sl accurate descrip- 
tion of the correlation hole of the extra particle is crucial 
for determining the band gap in crystals and sets thus 
the limitations of LDA-DFT. Naturally, various improve- 
ments have been suggested such as the o ptimized effec- 
tive potential method (see, e.g.,),'^2l2Sllll that allows for 
the use of orbital-dependent functionals in DFT. Orbital- 
dependent functionals were initially introduced within 
the exact exchange KS approach.''^ Alternatively, use can 
be made of Green's functions techniques, among which 
the most exploited is the GW approximation .'^ Yet, the 
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most enduring problem with the KS approach is that no 
systematic way has been developed to improve the func- 
tionals. 

Presently, advanced quantum chemical wave-function- 
based methods can address the electron correlation prob- 
lem at levels of high accuracy in small and medium-size 
molecular systems. These calculations work usually with 
virtual excitations of electrons out of molecular orbitals. 
Extending the application of wave-function-based meth- 
ods to solids is not straightforward since the delocalized 
nature of the Bloch orbitals would require the consid- 
eration of an infinite number of orbitals in order to in- 
corporate the relevant electron correlations. Needless to 
say, an ab initio calculation of the many-electron wave 
function based on Bloch orbitals scales quite unfavorably 
with the size of the system. For infinite systems we are 
concerned with, a relevant issue is to properly account for 
the translational symmetry while constructing the crys- 
tal many-electron wave function.'^'* ■^^ 

Despite the associated complexity, wave-function- 
based algorithms have been developed by several groups 
and applied to a series of semiconducting and insu- 
lating solid-state compounds. The electron correlation 
problem is addressed in these schemes either by second- 
order M0ller-Plesset theory (MP2)'* ^ ^ 38 39 40 coupled- 
electron-pair and coupled cluster methodsl ^ l ^ l ^ l 

By realizing that the correlation hole of an added elec- 
tron or hole is a local object, use can be made of lo- 
cal operators.l^^im The latter are most naturally asso- 
ciated with a set of real-space, localized Wannier or- 
bitals (WO's). Localized orbitals are usually obtained 
by applying various localization schemes to the canonical 
Bloch orbitals of t he self-con sistent-field (SCF) ground- 
state wave function.f^EiE^Esl Most often, su ch algor ithms 
are based on orbital localization procedurea^SHSIM! -^^jc];! 
were developed in order to facilitate correlation treat- 
ments in finite systems. The Foster-Boys localization 
criteriorP^ was successfully integrated in localization 
schemes for the canonical Bloch orbit als of per iodic sys- 
tems in several groups, see, e.g., Refs..'^I51Il2ll3l 

Localized orbital sets are extensively ex- 
ploite d in local correlat ion methods for 
solids! « l "' l " | i" | i^ | i4 | i^ l ^^ l ^« l ^" l ^'3 ^ A local approach 
to the computation of inter-atomic and intra-atomic 
correlations has been suggested by StoUhoff and Fulde.^^ 
Their formalism is based on a physically motivated 
selection of the relevant two-particle excitations by 
subdividing the orbital space into appropriate local 
orbital domains. This idea has been elaborated later on 
by Pulay and Saeb0, who developed correlated methods 
such as local M0ller-Plesset perturbation theory of 
second and higher (up to fourth) order without triple ex- 
citations and local coupled-electron-pair schemes. ^"^ ^'^ 
These local correlation approaches make use of a a 
priori division of the virtual orbital space into relevant 
subspaces, i.e., excitation domains. Although devel- 
oped initially for molecular systems, the local orbital 
schemes of Pulay and Saeb0 have been also adopted 



for crystalline compounds; see e.g..'i^Most recently, the 
local MP2 correlation approach of Pisani et al)^ has 
been proposed, which employs linear scaling techniques 
in order to facilitate the correlation calculation. This 
formalism was applied to the ground-state properties of 
various tetrahedral semiconductors. 

A local correlation method which can handle suc- 
cessfully excited-state properties of solids and goes be- 
yond the perturbatio nal schemes is th e local Hamilto- 
nian approach (LHA) ."i^IIlElIilllZIllIlIl The matrix ele- 
ments of the local Hamiltonian for the electron-addition 
and electron-removal quasiparticle states are computed 
in this approach by means of coupled-cluster or multi- 
configuration and multireference configuration interac- 
tion calculations^^ when the correlations are strong. The 
quasiparticle energy ban ds are then expr essed in terms 
of those matrix elements.E^ESElilliMlMl The computa- 
tional effort in the correlation calculation can be drasti- 
cally reduced by using the method of increments,'33T5?Illl 
designed initially to express the ground-state correlation 
energy as a fast convergent series of correlation contri- 
butions assigned to selected groups of localized orbitals. 
The incremental scheme has also become a valuable tool 
in constructing the scattering operator which is related to 
the ground-state wave functi on in a c umulant formulation 
of the many-body problem.'A^'^ZIll] The LHA has been 
applied to the determination of the correlation-induced 
corrections to the band structures of covalent semicon- 
ductors and insulators,*^ polymers,!^ model beryllium 
and hydrogen chains, and ionic insulators.!^ Extensive 
use of the incremental scheme was made for the accurate 
determination of the correlation contributions to the ex- 
cited state properties of diamond and silicon. SEHH! 

The starting point in our LHA is a ground-state peri- 
odic HF calculation. It provides a set of canonical Bloch 
orbitals which are then subject to an orbital localiza- 
tion procedure th at yields optimally localized Wannier 
orbitals; see, e.g.,.l^^*^*^^* ^ Alte rnative approaches were 
introduced by Shukla et aZ.EI^and Malrieu et al.^ who 
carried out the ground-state HF calculation in real rather 
than reciprocal space and obtained directly localized or- 
bitals. 

Ab initio many-body Green's function methods have 
also been de vised to com pute the energy bands of peri- 
odic systems.! "'^^ * ^" * ^^ * ^^ -*^ A local description of the corre- 
lation hole around an added electron or hole is vitilized in 
some of these formalisms too.^^ One of the key ad- 

vantages of such schemes is the fact that the self-energy 
operator is first evaluated exploiting local orbital sets and 
then transformed to a momentum representation. In ad- 
dition, the use of local orbitals prompts to employing the 
method of increme nts for the evaluation of the self-energy 
operator; see, e. g.,! 64 | 65| 

The key advantage of the wave-function-based meth- 
ods is the use of well defined and controllable approxima- 
tions; see, e.g., Ref. [56l They are amenable to systematic 
improvements that guarantee converged results for quan- 
tities such as binding energies and energy bands .'i^ 
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In the present study, we extend the application of the 
wave-function-based local Hamiltonian approach to cu- 
bic zinc-blende BN (c-BN). Due to its peculiar mechani- 
cal and physical properties, i.e., high bulk modulus, high 
thermal conductivity, and low dielectric constant, c-BN 
is a material of considerable technological interest In 
addition to experimental investigations of its electronic 
structure,®^ ^-^ a series of theoreti cal studies a re also 
available, b ased either on LDA-DFlEniniZlIMia or HF 
calculations.EmSl 

is a good candidate for extending 
previous studies on ionic o xidej^ and homonuclear co- 
valent compoundJi^EHIIIlIl the intermediate case of a 
partially ionic, heteropolar semiconductor. 

II. THEORY 

A. Formalism 

We begin with a brief review of the local Hamiltonian 
formalism and quasiparticle approximation. The starting 
point in our correlation treatment is a periodic HF calcu- 
lation that provides the self-consistent-field ground-state 
wave function |<J>gcF) and energy of the iV-electron sys- 
tem. Here, |<i>scF) is a single determinant wave function 
expressed in terms of the Hartree-Fock Bloch orbitals of 
the infinite system. The Bloch functions are eigenstates 
of the SCF part of the Hamiltonian, H — Hscf + -ffres, 
which is expanded in a Gaussian-type orbital (GTO) ba- 
sis, {/i(r)}, by using proper creation and annihilation 
operators aj^, a,J^^^^. 

^ = ^hija\^aj^+]^ 9i3kial„'^l^,ai„,aj„ (1) 

Here, hij and gijki are one-electron and two-electron ma- 
trix elements, respectively. Following the notations 
the compact index i of the set of operators {a]^} con- 
sists of an unit cell index / and an intracell index n. The 
partition of H into a SCF part, i?scFi and a residual in- 
teraction part, Hres, is essential within the framework 
of the local Hamiltonian approach and facilitates the ex- 
pression of the quasiparticle band energy as a sum of a 
SCF part and a correlation contribution. 

Next, we consider the HF electron-removal and 
electron-addition states, described in reciprocal k space 
as 

KJ.') - clj^sc.) and 
K'-^l') - Ck-..'l'i>scF), (2) 

where cj^^^^ creates an electron in the conduction-band 
Bloch state of momentum k, band index v, and spin 
CT, and Ck/i<T annihilates an electron in the valence-band 
Bloch state of momentum k', band index fi and spin a' . 
The essential idea of the LHA is to derive the correlation 
corrections to the HF band structure in terms of local 



matrix elements between explicitly constructed (iV-t-1) or 
(A^— 1) many-electron states. Due to the local character of 
the correlation hole, these Hamiltonian matrix elements 
are obtained in practice from calculations on sufficiently 
large finite clusters. To accomplish this local description, 
the one-electron Bloch orbitals are expressed as a super- 
position of Wannier orbitals associated with the canon- 
ical HF valence and conduction bands. The operators 
•^ki/cr are expanded in terms of operators for the Wannier 
orbitals wL as 

^L.-^J2^-^-(^)<n.e'''-''\ (3) 
^ ° n,Ri 

where Ri is the lattice vector of the unit cell / at which 
the Wannier function |w„o-(R.i)) is centered. A^o is the 
number of unit cells. An analogous expression also ap- 
plies for Ck'^o-'- Hence, becomes 

= 7^ E a.«(k)e^'^-^^<„J<i>scF) (4) 

and a similar expression holds for |*&^^^/). To this end, 

l*^ki«T^) ^'^^ subject to a transformation, 

analogous to the Wannier transformation, which yields 
local electron- addition and electron-removal one-particle 
configurations 

I'J'r;™.') = W^.ma'\^scr). (5) 

The energies of the conduction bands are expressed in 
terms of Hamiltonian ma trix elemen ts between those lo- 
cal configurations l^^+i^^IlMIlH 

4^". = 5]^a.„(k)a:„,(k)e'''-^'i/|,-„„„ (6) 

R.I nn' 

where 

Hr^-r^^uu' = (*scFk'RK«'ai?w;J^,„<,|«'scF) 

E^- (7) 

These matrix elements are uniquely defined. Here, Eq'-'^ 
is the Hartree-Fock ground-state energy. An expression 
analogous to Eq. ^ applies for the energies of the va- 
lence bands. For the valence-band states, the matrix el- 
ements have the form 

^Rj-RL,mm' " " RjRl "mm' -^^0 

- (*SCF|wil^,„V'^WR,,„,.|^sCF). (8) 

The local Hamiltonian matrix elements in Eqs. ([t]) 
and (Isl are ext racted from calculations on large enough 
clustCTS.™Ill! 

Next, we consider the inclusion of electron correla- 
tions in the band structure calculations. For relatively 
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weakly correlated electron systems, such as boron nitride, 
we may adopt the quasiparticle approximation.^* 
Within this approximation, the life time of the excita- 
tions is neglected. Satellite structures are also beyond 
the quasiparticle picture. In analogy with the HF ap- 
proximation, the correlated conduction and valence en- 
ergy bands are expressed as 



Eq and 



(9) 



where Eq is the energy of the correlated ground state 



and 1^1' 



and l^'^^J/) are the correlated counter- 



parts of |$[(.^^ ) and |<i>i(.7icr')' respectively. If a Wannier- 
likc transformation is applied to the wave functions 



N+l\ 



and l^-" 



tions 



N+l 



and l^^-i 



the locally correlated wave func- 
are obtained. The latter 



are viewed as derived from the correlated ground-state 
wave function, after an electron-addition to the conduc- 
tion band or an electron-removal process from the valence 
band has taken place and the subsequent response of the 
crystalline surroundi ngs to the add itional charge is ex- 
plicitly accounted for.™ElIiIllZllII 

In analogy with the HF case, the energy bands are ex- 
pressed in terms of real-space matrix elements between 
localized (-/V±l) states. In this case, however, the matrix 
elements incorporate electron correlation effects. The 
conduction energy bands are given by 



cor7' 
<^ki/(T • 



(10) 



The expression for the correlated valence bands is anal- 
ogous. The expressions for the energy bands can also 
be written in terms of cumulants, which insures size- 
extensivity of the different correlation energy contribu- 
tions. The cumulant formulation has been discussed in 
detail in Refs. (THlandllTl 

In real space, we can represent the correlated (A^+l) 
electron states as 



I* 



AT + l 

R.Tn(T 



e^|$ 



Rino-/ ' 



(11) 



and similarly for the localized and correlated {N—Vj elec- 
tron states. The operator S is constru cted from a selected 
set of excitation operators; see, e. g.j.EMZIilllZI Por rela- 
tively weakly correlated systems, a reasonable ansatz for 
the opera tor S inclu des one- and two-particle excitation 
operators,™l23111 



S 



/ ^ / ^ / ^ '«.jK,kJ^L^IVI 



where the indices m, m' refer to the occupied Wannier 
orbitals of the (-/V+l) electron system and a, 6 denote 
virtual orbitals centered at Rj, Rk- The expression for 
the operator S for the (-/V— 1) electron system is analogous. 

The operator S has a twofold function. The first 
term in Eq. ( 12 1 generates the relaxation and polar- 



ization cloud around an extra electron, created by the 
composite operator w\^^^^^'WYl^nac\^^^ in the conduction- 
band Wannier orbital |w„ct(Ri))- Analogously, a relax- 
ation and polarization cloud is formed around a hole cre- 



,j mo-' U^RjniCT' "^k'/jcr' 



in the valence-band WO 



ated by 

|fi;„iCT'(Rj))- The second term in Eq. \12\ accou nts for 
the so-called loss of ground- state correlation! ^ ''I^^ l We de- 
scribe first the procedure for the computation of the first 
term. 

On-site and short-range inter-atomic relaxation and 
polarization effects are computed by means of separate 
SCF orbital optimizations for a large but finite region 
around the a dditio nal charge. The frozen local hole 
approximatioiJSEIl is applied, i.e., the Wannier orbital 
where the extra hole or electron resides is kept frozen in 
the SCF calculatioii2ai for the (iV±l) excited state. This 
approximation is well suited if the electrons are mod- 
erately correlated. It has been demonstratecf^ that up 
to first-order perturbation theory, the SCF orbital re- 
laxations are equivalent to the one-particle excitations 
around the frozen Wannier orbital, see the first term in 



Eq. (12). The corresponding correlated wave functions 



Rjmcr'/ ^^"^ I^Rino-/'^^^'^^ ^^'^^ identical up to first or- 



der to the optimized SCF states |*&Rj„cr'/ 



and |<i> 



N+l 

Ri71(t/ 



respectively. Note, that the creation and annihilation op- 
erators in Eq. (12 1 refer to spin- orbitals. In practice, a 



spatial orbital set is exploited and therefore we also allow 
for the relaxation of the orbital accommodating the extra 
electron or hole within the computed relaxed surround- 
ing in order to account for the spin degree of freedom; see 
Section IV. Wave functions which are expressed in terms 
of relaxed orbitals within the nearby surroundings and 
readjusted hole or added-electron orbitals are denoted in 



the following as \^'^^^,,) and |$r,„,,. 

The long-range polarization of the crystal brings also 
large corrections to the diagonal matrix elements of the 
local Hamiltonian and is estima ted within the approxi- 
mation of a dielectric continuum.iSlilllZllIl 

Next, we consider the loss of ground-state correlation. 
This effect is related to correlation contributions that are 
absent in the {N ± 1) electron system but present in the 
A^-electron ground state. This has to do with the fact 
that some excitations involving the removed electron or 
the orbital where the extra electron is placed are blocked 
in the {N ± 1) electron system. 

To account for electron correlation effects beyond re- 
laxation and polarization, an important part of which 
is the loss of ground-state correlations, multi-reference 
(MR) single and double configuration interaction (SDCI) 
calculation^^2lin! are carried out for each 1$ 



(12) i\^^:L')) and ||.^tV) (1$^+!,)) state. The CI wave 



' R.J mcT' 
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functions, and l^^'^^^), are constructed within 

an orbital space consisting of the relaxed doubly occu- 
pied SCF orbitals, the singly occupied, [^^^'(Rj)) or 
|w„(j(Ri)), and the set of virtual orbitals. Within the 
quasiparticle approximation, the CI expansions contain 
only configurations in which the hole in |i(;„i(j'(Rj)) or 
the extra electron in |w„(j(Ri)) is never promoted to an- 
other orbital. Differential correlation effects arising from 
having a different number of electrons in the N and {NzLT) 
configurations are also assessed by means of CI calcula- 
tions. 



B. Calculation of matrix elements 

The required matrix elements can be obtained from 
cluster calculations, as pointed out above. The crys- 
talline environment around such a cluster is described by 
an embedding potential which is extracted from the pe- 
riodic Hartree-Fock calculation. This embedding scheme 
was proposed by Birkenheuer et al., see RefslT4l andSl] 
and requires a prior localization of the canonical HF or- 
bitals of the A^-electron periodic system. The practical 
evaluation of the embedding potential is facilitated by 
the knowledge of the crystal Fock operator Fcryst and the 
cluster Fock operator Fciust[Pciust\- The latter is related 
to the density operator Pdust = Ylma"^^"^ \wmcr){wmcr\, 
where \wm(j) are the localized occupied orbitals within 
the cluster. 

The cluster is subdivided into an active region and a 
buffer region. The occupied and low-lying conduction- 
band orbitals of the active region are explicitly correlated 
in the post-HF calculation, while those of the buffer re- 
gion are always kept frozen. The role of the atoms in the 
buffer region is to insure an accurate representation of the 
longer-range tails of the Wannier orbitals, |w,„(j'(R-j)) 
and |w„cr(Ri)), centered at the atomic sites J and I of 
the active region; se^SEU and also Section IV. 

The orbital set associated with the finite cluster is 
generated from the original crystal Wannier orbitals. 
The latter are o btained via a Wannier-Boys localiza- 
tion schemd^mH for the core, valence, and low-lying 
conduction-band states. Further, the crystal WO's are 
projected onto the s et of atomic basis functions assigned 
to the cluster sites.'SEH These projected Wannier or- 
bitals, (Rj)) and |it;^^(Ri)), do not constitute an 
orthonormal orbital set. Therefore, a series of groupwise 
orthonormalizations are carried out for the core, valence, 
and low-lying conduction-band orbital groups, both for 
the active and buffer regions. The sets of orthonormal- 
ized occupied and low-lying conduction-band orbitals are 
denoted as |wJ^^/(Rj)) and |w"^(Ri)). 

Additionally, a set of virtual orbitals is needed in or- 
der to construct the complete variational orbital space 
for the correlation calculation. To generate those vir- 
tual orbitals, a modified version of the projected atomic 
orbital (PAO) formalism of Hampel and Wernei'^ is 
adopted.^'' The idea for using projected atomic orbitals 



as a local virtual basis was initially proposed by Pulay.l^ 
The virtual PAO's are obtained from the atomic or- 
bital basis of the active cluster region after projecting 
out the contribution of the orthonormalized, occupied 
and low-lying conduction-band orbitals |w^^/(Rj)) and 
|fi;^^(Ri)) via a Schmidt orthogonalization scheme; see 
also Ref.TSl The PAO's are subsequently Lowdin or- 
thonormalized. The final variational orbital space of 
the embedded cluster consists of the core, valence-band, 
and low-lying conduction-band orbitals |i(j''jg.,(Rj)) and 
|wJ^„(Ri)) together with the set of virtual PAO's. The 
local electron-addition and electron-removal configura- 
tions, \'^^na) ^'^'^ l*^R7mo-')' ^'^^ defined within this vari- 
ational orbital space. 



III. CRYSTAL STRUCTURE OF C-BN AND 
COMPUTATIONAL INFORMATION 

Boron nitride belongs to the group of III- V binary com- 
pounds and exhibits three solid state phases. At room 
temperature and normal atmospheric pressure it crys- 
tallizes in a hexagonal structure, which transforms to 
dens er pac ked zinc-blende or wurtzite structures under 
stati(P^Ell or dynamic compression respectively. The 
cubic zinc-blende form is face-centered (fee), with the 
FAim space group symmetry. The corresponding lattice 
constant is 3.615 A .^^ The local point-group symmetry 
for both B and N is T,;, with a four site nearest-neighbor 
coordination. 

GTO basis sets of triple-zeta quality augmented with 
polarization functions were applied for the N atoms. We 
used as a starting point the N 7-311G GTO basis set de- 
rived by Dovesi et al^^ for crystalline LiN. We added to 
Dovesi's basis set a single-Gaussian d polarization func- 
tion with an exponent of 0.8^^ ^'^ and reoptimized the 
exponents of the outermost two sp shells. The exponents 
of the reoptimized outermost sp shells are 0.4523 and 
0.220241. 

Two different GTO basis sets were used for the B 
atoms. The first is a triple-zeta B 6-31 IG* basis set which 
is derived from the B 6-21G* GTO basis set constructed 
for BN by Orlando et aZ.. ESI This 6-21G* basis set was 
suitably modified for the purpose of obtaining a triple- 
zeta basis set. The optimized values of the exponents 
and contraction coefficients of the valence 2sp shell and 
the exponents of the polarization single-Gaussian isp and 
4sp shells are listed in Table |V] of the Appendix. 

The second basis set is a double-zeta 6-41G basis set 
which is designed starting from the 6-31 IG* basis set 
given in Table |Vj This smaller basis set is employed for 
calculations on very large clusters, for which the use of 
the 6-311G* basis is at present computationally unfeasi- 
ble. Those large clusters are selected in order to compute 
additional increments. As discussed in the following sec- 
tion, the associated basis set effect is either negligible or 
very small. To obtain the 6-41 G basis set, the d polar- 
ization function is removed and the single-Gaussian 3sp 
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function is contracted with the three primitives of the 2sp 
shell. The exponents and contraction coefficients of the 
2sp shell and the exponent of the new single-Gaussian 
3sp function are reoptimized. Their values are provided 



in Table VI of the Appendix. 

The periodic HF calculations are performed with the 
CRYSTAL program package .1221 crystal provides in addi- 
tion a Wannier-Boys orbital localization module for gen- 
erating real-space, localized Wannier functions.!^ The set 
of projected Wannier orbitals for the finite cluster and the 
matrix representation of the crystal Fock operator Fcryst 
in terms of the cluster atomic basis functions are both ob- 
tained with the CRYSTAL-MOLPRO interface program.!^ 
The embedding potential for the finite cluster is gen- 
erated with MOLPRcl^ by substracting from Fcryst the 
cluster Fock operator Fdust [Pciust] , the latter being com- 
puted from the density matrix associated with the set of 
occupied, projected WO's. All subsequent calculations 
for studying the effects of electron correlations are per- 
formed with the MOLPRO package.!^ 



IV. CORRELATION INDUCED CORRECTIONS 
TO THE BAND STRUCTURE OF ZINC-BLENDE 

C-BN 

The HF energy bands of c-BN, obtained with the 
triple-zeta basis sets described in the previous section, 
are plotted in Fig. [T] The HF energy per unit cell is 



1.00 



0.50 



3 
LU 



.0.00 



-0.50 



-1.00 




W K 



W 



FIG. 1: Hartree-Fock band structure of c-BN. The core N Is 
and B Is bands are not shown in the figure. 



—79.268 a.u., close to the value obtained by Euwema et 
ai, —79.250 a.u..-^ The upper valence bands of c-BN 
have N 2p character, with small contributions from the 
B 2s and 2p states. The low-lying conduction bands have 
predominant B 2s and B 2p character, with sizeable ad- 
mixture from the N 2s, 3s and 2p,3p functions. The fun- 



damental band gap of the system is indirect, F 



15 



At the Hartree-Fock level this quantity is strongly over- 
estimated, i.e., 13.62 eV, and therefore 7 eV larger than 
experimental values dedu ced fr om soft x-ray experiments 
at the B and N if-edgesPEl and optical absorption^^l 
measurements. Those values are 6.0±0.5 and 6.4±0.5 
eV, respectively. The energy separation between the N 
2p and lower-lying N 2s valence bands is about 6 eV. 




FIG. 2: One of the four conduction-band sp'^-fike Wannier 
orbitals with predominant B 2s, 2p character after projection 
onto a [B13N28] cluster. 



Wannicr-Boys transformations are carried out sepa- 
rately for the core, valence, and low-lying conduction 
bands. The lowest four conduction bands of BN are sepa- 
rated by a small gap from the higher virtual Bloch states. 
Due to this finite energy gap, no numerical problems oc- 
cur in the Wannier-Boys transformation for the set of 
low-lying conduction-band states. The resulting WO's 
associated with these lowest conduction bands turn out 
to be a set of four sp^-like hybrids, each oriented along 
one of the B-N segments of a BN4 tetrahedron. They 
have large weight at the nearest-neighbor N sites and 
strong antibonding character with respect to the nearest- 
neighbor N 2s and N 2p orbitals. Such a Wannier orbital, 
projected onto a finite cluster, is plotted in Fig. |2] 

In a next step, the projected WO's are subject to a 
Pipek-Mezey transformation.!^ In simple cases, this pro- 
cedure transforms a set of sp hybrids into Cartesian type, 
s and p functions, as found for example in the case of the 
conduction-band Wannier orbitals of MgO.*^ Here, the 
three degenerate orbitals, each oriented along one of the 
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C2 axes, and the fourth, lower-energy WO obtained after 
the Pipek-Mezey transformation bear little resemblance 
to Cartesian functions. Each p-like lobe of the three de- 
generate orbitals at the B site of a BN4 tetrahedron is 
antibonding with respect to the two N ions adjacent to 
that lobe, but has sub stan tial, U-shaped bonding tails 
at the other two N site^^^H, xhe fourth function is also 
strongly deformed as compared to a pure s function. In 
the immediate neighborhood of a B site, it is tetrahe- 
drally shaped. These projected WO's, obtained through 
the Pipek-Mezey transformation, are plotted in Figs. [3] 
and m 




FIG. 3: B 2p-like conduction-band Wannier orbital after pro- 
jection onto a [B13N28] cluster. 



In contrast to the conduction-band Wannier functions, 
the Wannier orbitals associated with the nitrogen 2s and 
2p valence bands are well localized around the N sites, 
with very small tails at the nearest boron neighbors. The 
Pipek-Mezey transformation yields a set of somewhat de- 
formed, Cartesian-like 2s, 2px, 2py, and 2pz functions. 
One of the 2p-like components and the deformed N 2s 
orbital are plotted in Fig. |5] and Fig. [6) respectively. 

Worthwhile to mention, the norms of the projected 
valence-band WO's, centered within the active regions 
of the embedded clusters, are always larger than 0.99 of 
the initial HF WO's obtained with crystal. The corre- 
sponding norms for the projected conduction-band WO's 
are always larger than 0.94. 



A. Diagonal matrix elements of the local 
Hamiltonian 

1. Short-range relaxation and polarization effects 

We start by considering the short-range correlation- 
induced corrections to the diagonal matrix elements of 




FIG. 4: B 2s-like conduction-band Wannier orbital after pro- 
jection onto a [B13N28] cluster. 




FIG. 5: N 2p-like valence-band Wannier orbital, projected 
onto a [N13B28] cluster. 



the Hamiltonian in Eqs. ^ and ([s]). We construct 
first the correlation cloud of a hole introduced in the 
valence bands. To evaluate the short-range part of this 
correlation cloud, i.e., the on-site relaxation as well as 
the relaxation and polarization effects associated with 
the first (nearest-neighbor, nn) and second (next-nearest- 
neighbor, nnn) nitrogen coordination shells around the 
N site where the hole resides, we designed three differ- 
ent clusters: [N55B28], [N71B36], and [N63B44]. The first 
nitrogen coordination shell contains 12 sites whereas the 
second shell consists of 6 N atoms. 
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TABLE I: Correlation-induced corrections (in eV) to the diag- 
onal Hamiltonian matrix elements for the valence-band N 2s 
and 2p hole states. Negative corrections indicate an upward 
shift of the valence bands. 



A H™(0) 




N 2p 


N 2s 


On-site orb. relaxation 


-1.29 


-1.64 


nn(B) orb. relaxation 


-0.01 


-0.01 


nn(N) orb. relaxation 


-1.57 


-1.63 


nnn(N) orb. relaxation 


-0.14 


-0.14 


Hole orb. relaxation 


-0.18 


-0.12 


Loss of ground-state corr. 


1.05 


1.08 


Long-range polarization 


-1.53 


-1.53 


Total corr. correction 


-3.67 


-3.99 



FIG. 6: N 2s-like valence-band Wannier orbital, projected 
onto a [N13B28] cluster. Due to the local environment, it is 
distorted to a tetrahedral-like shape. 




FIG. 7; Illustration of the [N71B36] embedded cluster. The 
active region of the cluster consists of a [N17B4] fragment. 
The N sites in that region are depicted as black spheres of 
different size, depending on their distance from the central N 
site. The latter is shown as the largest black sphere. N atoms 
in the buffer zone are shown as small, grey spheres. The four 
active B atoms are drawn as large white spheres, whereas the 
buffer B sites are shown as small white spheres. 



The [N55B28] cluster is employed for computing the 
on-site and nn relaxation and polarization effects. The 
active region, [N13B4], of the [N55B28] cluster includes a 
central N site, sketched as the largest black sphere in Fig. 
|7] the four nn B atoms, depicted as large white spheres, 
and the first shell of 12 nn N atoms. These 12 N nn's are 
drawn as average-size black spheres in Fig. [7] The buffer 



region around the active subunit contains all B and N 
neighbors in the first two coordination shells of the 12 
outer N sites of the active region. Such a buffer region is 
large enough to ensure a good representation of the long- 
range tails of the Wannier orbitals centered in the active 
subunit. In this set of calculations, we exploit the 6-41G 
basis set for the B atoms (see Section III) because the use 
of the 6-311G* basis set is computationally unfeasible for 
such a large cluster. The nitrogens are described by the 
7-311G* basis set. 

Separate restricted open-shell HF (ROHF) calculations 
are carried out for the N 2s and N 2p {N — 1) electron 
states. First, the local electron-removal states I'&R^mcr') 
are constructed according to Koopmans' theorem. In a 
second step, the hole is placed in one of the N 2s or 2p 
WO's at the central N site of the [N55B28] cluster and 
the other valence orbitals within the active region are 
allowed to relax and polarize in response to the electron- 
removal pr ocess. The Wannier orbital with the hole is 
kept frozen.f2Ii^ We construct thus l^^J^a') *s a first 
approximation to the correlated wave function l^'ajmcr'); 
see also Section II. 

The corrections due to on-site orbital relaxation are 
found to be substantial, 1.3-1.6 eV for the N 2s and 2p 
hole states (first line in Table |l|. The effect is stronger 
for the N 2s hole state because in this case the on-site 
relaxation involves all six electrons in the highly polariz- 
able N 2p shell. The relaxation of the nn boron Is core 
shells is only a small fraction of 0.01 eV of the overall 
correlation-induced shift of the valence bands, as shown 
on the second line in Table IT] Therefore this effect is not 
discussed further. The relaxation and polarization of the 
nn N orbitals bring corrections of similar magnitude for 
both types of hole states, about 1.6 eV (third entry in 
Table |l|. The relaxation of the nitrogen Is core shell is 
vanishingly small. For this reason, all N core shells are 
kept frozen throughout the calculations. 

Some comments regarding the basis set used for the 
B atoms are here in place. Basis set effects were investi- 
gated by additional calculations for the on-site relaxation 
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and polarization using a [N13B28] cluster. The active 
region, [NB4], of the [N13B28] cluster contains only the 
central N site and the four nn B atoms. The first shell 
of twelve nn nitrogens around the central N site is now 
part of the buffer region. ROHF calculations are per- 
formed for the N 2s and N 2p (N—t) electron states of 
the [N13B28] cluster using for the borons either 6-41G or 
6-311G* basis sets. The basis set effect for the on-site 
orbital relaxation is found to be negligibly small, <1%, 
for both the N 2s and and 2p hole states. 

In a next step, the contribution of the nnn nitrogens 
to the overall relaxation and polarization has been in- 
vestigated. Since wave-function-based calculations for a 
cluster containing all nn and nnn N sites are computa- 
tionally unfeasible even with the double-zeta basis set, 
we designed clusters where only part of the six nnn N 
sites were included in the active region. For one of these 
clusters, [N71B36], the active region [N17B4] includes the 
active region [N13B4] of the [N55B28] cluster and four ni- 
trogens from the second coordination shell of the central 
N site. These four nnn N's lie along the and 
axes. Analogously, the active region [N15B8] of the other 
selected cluster, [N63B44], contains the [N13B4] fragment, 
the two nnn N atoms along the C2X axis plus the four 
borons bridging those two nnn nitrogens with the [N13B4] 
subunit. The calculations are carried out using the 6-41 G 
basis set for B and the 7-311G basis for N^^"*. 

The overall correlation correction associated with or- 
bital relaxation and polarization effects at the six nnn ni- 
trogens is listed on the fourth line of Table |lj It amounts 
to -0.14 eV and was obtained as a sum of the correlation 
corrections originating from the orbital relaxation and 
polarization at the four nnn N in [N7iB3g] and the two 
nnn N in [N63B44] . To confirm that those corrections are 
indeed additive, we performed additional calculations on 
[N7iB3g], in which we allowed for the orbitals of only two 
of the four nnn N sites to relax and polarize. The relax- 
ation effect for these two N atoms amounts to approxi- 
mately -0.045 eV, which represents indeed one third of 
that -0.14 eV mentioned above. It is interesting to note 
that the correlation corrections associated with the nnn 
N sites are independent of the angular momentum of the 
orbital where the hole resides. This finding indicates that 
the interaction between the extra charge, placed at the 
central N, and the induced dipole at those nnn N sites is 
a monopole-dipole type interaction. 

At last, we calculated the correction due to the relax- 
ation of the orbital where the hole resides (fifth line in 
Table |l]) . We accounted in this manner for the spin degree 
of freedom after the electron-removal. The hole orbital 
localized at the central N is now allowed to relax in the 
presence of nearby valence orbitals, on-site, nn, and nnn, 
that were allowed to relax in the first set of calculations. 
The nearby nitrogen orbitals around the hole are kept 
frozen during this optimization of the singly occupied or- 
bital. 

Worthwhile to mention, calculations of the on-site or- 
bital relaxation and polarization correction, AHJ'^^^(O), 



TABLE IL Correlation-induced corrections (in eV) to the di- 
agonal Hamiltonian matrix elements for the conduction-band 
B 2s and B 2p electron-addition states. Negative corrections 
indicate shifts to lower energies of the conduction bands. 

A H„„(0) 

B 2s-like B 2p-like 



nn(N) orb. relaxation -1.03 -1.14 

nnn(N) orb. relaxation -0.42 -0.45 

Added-elec. orb. relaxation -0.14 -0.18 

Loss of ground-state corr. 0.60 0.50 

Long-range polarization -1.77 -1.77 

Total corr. correction -2.76 -3.04 



exploiting either [N13B28], [N55B28] or [N63B44] clusters 
indicate that t his quantity is converged with respect to 
the cluster siz^l^SJ. Similar observations hold for the re- 
laxation and polarization corrections associated with the 
nn N sites. 

Having computed the relaxation and short-range po- 
larization effects on the diagonal matrix elements for the 
(A^— 1) valence-band states, we discuss next the low energy 
(A + 1) conduction-band states. These electron-addition 
states imply configurations with an extra electron into 
the B 2s-like or B 2p-like conduction-band WO's. We 
employed a [B31N68] cluster to investigate such configu- 
rations. The active region [BNig] of the [B31N68] cluster 
consists of a central B site and the first two nitrogen 
coordination shells, containing four and twelve N atoms, 
respectively. The buffer region includes again all B and N 
sites within the first two coordination shells of the outer 
12 active N atoms. We used in these calculations the 
6-41G basis set for B and the 7-31 IG basis set for N. 

In analogy to the (A^— 1) valence-band states, we have 
carried out ROHF calculations for the lowest-energy (A-|- 
1) states keeping the singly-occupied B 2s-like or B 2p- 
like orbitals frozen.l^ The valence 2s and 2p orbitals of 
the N atoms within the active region [BNig] are allowed 
to relax and polarize in response to the electron-addition 
process. 

The on-site relaxation effect associated with the B Is 
core orbital is found to be vanishingly small. The relax- 
ation of the nn nitrogen atoms gives rise to a downward 
shift of 1.0-1.1 eV for the B 2s and 2p conduction bands, 
see Table |TT] Relaxation effects at the nnn N sites bring 
an additional downward shift of about 0.45 eV of the 
center of gravity of the four low- lying conduction bands. 

In a next step, the B 2s-like or B 2p-like orbital with 
the extra electron is also allowed to relax, while keep- 
ing the other valence and low-lying conduction-band or- 
bitals frozen. The reason for freezing the latter, i.e., the 
unoccupied B 2s and 2p-like components, is to prevent 
delocalization of the extra electron. For the same rea- 
son, in the case of the B 2p added-electron states, the 
next four higher-energy virtual orbitals are kept frozen 
as well. Those are linear combinations of s and p virtual 
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orbitals centered at the nearest N and B sites. Worth- 
while to note is the similar magnitude of this correction 
as compared to the relaxation correction found for the N 
2s, 2p hole orbitals (see Table . 

An analysis of basis set effects similar to that per- 
formed for the {N— 1) states was also carried out for the 
electron-addition states. Additional calculations of the 
relaxation effects associated with the nn N atoms were 
performed on a [B13N28] cluster. The active region of 
this cluster contains the central B site and the four nn N 
sites. The second nitrogen coordination shell of the cen- 
tral B (i.e., nnn N atoms) is now included in the buffer 
region. ROHF wave functions were constructed for the B 
2s and 2p (A^-l-l) states using for the boron and nitrogen 
atoms either B 6-41G and N 7-31 IG or B 6-31 IG* and 
N 7-311G* basis sets. As compared to the {N—1) states, 
the effect is somewhat larger for the (A^-l-l) states, 6 to 
7 %. The correlation-induced correction associated with 
the short-range relaxation at the four nn N sites (first line 
in Table includes this contribution of 6-7 % related to 
basis set effects. 

To summarize this subsection, the on-site orbital relax- 
ation and short-range relaxation and polarization effects 
at the nn and nnn nitrogen sites in response to the extra 
charge, hole or electron, give rise to corrections of several 
eV to the HF band gap. 

2. Long-range polarization effects 

A substantial contribution to the reduction of the HF 
gap also arises from long-range polarization effects. A 
simple estimate of the long-range polarization contribu- 
tion to the shift of the energy bands is assessable through 
the computation of the classical polarization energy of a 
dielectric medium outside a sphere of radius R, beyond 
which the dielectric response of the crystal reaches its 
asymptotic value, eg^'': 



This approximation is applicable because the interaction 
of the extra charge with the surroundings beyond a given 
R has predominantly electrostatic character. The re- 
laxation and short-range polarization effects within the 
sphere were already computed at the ab initio level and 
are summarized in Tables U and |lll 

Using the experimental static dielectric constant, 
£0=6.8,1^ we find corrections of few eV to the diagonal 
matrix elements for the valence- and conduction-band 
states. Those values are also included in Tables HI and ITTl 
The cuttoff radii R for the {N-1) and (N+l) excited 
states, are each obtained as the average of the radii of 
the second and third nitrogen coordination shells around 
the central N or B sites. The corresponding values are 
4.021 A for the valence-band hole states and 3.468 A for 
the conduction-band electron-addition states. The long- 
range polarization corrections ^H'^n7n(nn)(^) then - 



1.53 eV for the N 2s, 2p hole states and -1.77 eV for the 
B 2s, 2p electron-addition states. 

It is instructive to compare the ab initio result for the 
polarization effect at the nnn N sites with that obtained 
within the dielectric continuum approximation. We re- 
call that around each N site there are 12 N atoms in the 
first nitrogen coordination shell (nn N) and 6 N sites in 
the second N coordination shell (nnn N). The ab initio 
value for the contribution of the nnn N atoms to the over- 
all correlation-induced correction to the diagonal matrix 
elements for the N 2s, 2p hole states was found to be 
-0.14 eV, see Table |T] and the paragraphs above. A sep- 
arate estimate for this quantity can be obtained within 
the dielectric continuum approximation. In this approx- 
imation, the correction to the diagonal matrix elements 
due to polarization effects at the second shell of N neigh- 
bors is 5E2 = S2^e^{^ - with i?i=3.086 A, 

i?2=4.021 A , and S2 =0.47. Here, S2 is a scaling factor 
which represents the ratio between the density of N sites 
in the spherical shell enclosed between Ri and R2 and the 
density of N atoms in the unit cell, see also.^i^ The cor- 
rection 6E2 is found to be -0.21 eV, which is about 30% 
larger than the ab initio value of -0.14 eV. This difference 
between the ab initio estimate and the value deduced on 
the basis of the continuum dielectric model is not very 
surprising. A better agreement between the two values 
would imply that a continuum approximation may be al- 
ready made for the polarization beyond the first shell of 
N neighbors. This can not be expected to hold because 
at such short distances the associated charge distribution 
and polarization is not that of a continuum. Therefore, 
the errors related to the continuum dielectric approxima- 
tion axe larger for small R's. 

3. Loss of ground-state correlation 

An important contribution to the correlation-induced 
corrections to the diagonal matrix elements also arises 
from differential correlation effects, i.e., loss of ground- 
state correlation. To study such effects, we construct 
the correlated wave functions I^r^wkt') I^RitiCT) 
means of SDCI calculations. 

The localized character of the valence-band WO's al- 
lows for the use of a [N13B28] cluster, with a single N 
atom in the active region [NB4], for studying correlation 
effects beyond charge relaxation and polarization. The 
basis sets exploited in the SDCI calculations are 6-311G* 
for the B and 7-311G* for the N atoms. 

The SDCI wave functions for the N and (TV — 1) elec- 
tron states are constructed by correlating the 2s and 2p 
orbitals of the central N atom in the [NB4] kernel. For 
the (iV— 1) electron configuration, the reference wave func- 
tions for the different N 2s and N 2p hole states l^^J^a') 
are each expressed in terms of individually optimized or- 
bital sets. The orbital with the hole is also relaxed, as 
described above. In the framework of the quasiparticle 
approximation, the SDCI wave function for the (iV — 1) 
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state should be expanded in terms of configuration state 
functions (CSF's) for which the occupation of the N 2s 
or 2p hole spin-orbital is kept frozen. This is referred to 
as the frozen local hole approximation.^ In practice, we 
work with spatial orbital rather than spin-orbital sets. In 
MOLPRO, for example, the configuration selection is done 
for spatial orbitals and hence the SDCI wave function, 
constructed by correlating the singly occupied N 2s or 
2p orbital with the doubly occupied N 2s, 2p orbitals, 
may contain CSF's for which the occupation of the hole 
orbital is no longer maintained. The contributions of 
such CSF's to the total correlation energy for the (-/V— 1) 
state account for efi^ects beyond the loss of ground-state 
correlation and the quasiparticle picture. The associated 
excitations account for satellite structures in the pho- 
toionization spectra. To assess the various correlation 
contributions, in particular, the loss of ground-state cor- 
relation, we analyzed the correlation energy for the (A^— 1) 
states in terms of intra- and inter-pair contributions of 
the N 2s and 2p orbitals. 

Consider one of the three N 2p hole states. In a first 
step, we included in the correlation treatment all sin- 
glet and triplet electron pairs excited from the doubly 
occupied N 2s and 2p orbitals and the singly occupied 
N 2p hole orbital. The associated SDCI correlation en- 
ergy of about -4.00 eV is partitioned into contributions 
of excited internal, singly external, and doubly external 
configurations. The terminology used for the different 
types of configurations follows that adopted by Werner 
and Knowles in Ref. i79i The excited singly external con- 
figurations are related to one-particle and semi-external 
two-particle excitations for which one electron is excited 
from a singly or doubly occupied orbital to a virtual (ex- 
ternal) orbital. The excited doubly external configura- 
tions result from two-particle excitations for which both 
electrons are promoted to the virtual orbital space, i.e., 
external two-particle excitations. 

We analyzed the different contributions of these ex- 
cited internal, singly external, and doubly external con- 
figurations to the total SDCI correlation energy. With 
the present variational orbital space, excited doubly ex- 
ternal configurations do not change the occupation of the 
hole orbital. This means that the intra- and inter-pair 
correlation contributions of the N 2s and 2p orbitals orig- 
inating from such configurations are a relevant part in 
the computation of the loss of ground-state correlation. 
Contrary to the doubly external configurations, the ex- 
cited internal and the singly external configurations in 
which the occupation of the hole orbital is altered ac- 
count for effects beyond the quasiparticle picture and 
frozen local hole approximation. Hence, the associated 
correlation contributions should not be considered in the 
evaluation of the loss of ground-state correlation. Our 
calculations indicate that with the present choice of the 
cluster, the total SDCI correlation energy of -4.00 eV 
contains a negligibly small contribution from excited in- 
ternal configurations, less than -0.001 eV. The overall 
correlation contribution arising from the excited singly 



external configurations is -0.53 eV. We are interested in 
determining the fraction due to configurations in which 
the occupation of the hole orbital is preserved. An anal- 
ysis of the SDCI wave function indicate comparable co- 
efficients for the excited singly external CSF's in which 
the N 2p hole orbital is either doubly or singly occupied. 
The configurations with a doubly occupied N 2p "hole" 
orbital arise from so-called semi-external two-particle ex- 
citations. Since the composition of the SDCI wave func- 
tion is not sufficient to determine the relevant part of the 
net effect, we carried out further the following analysis. 

In a second step, we computed all various intra- and 
inter-pair contributions of the N 2s and 2p orbitals by 
correlating in separate SDCI calculations different or- 
bitals and orbital pairs, respectively. We constructed 
thus SDCI wave functions where either one orbital, a 
pair of orbitals, or a combination of any three orbitals 
from the N 2s, 2p orbital set were correlated. These 
wave functions were carefully examined to determine the 
singly external CSF's with large coefficients in the CI ex- 
pansions. We found that the correlation energy due to 
excited singly external configurations in which the occu- 
pation of the N 2p hole orbital is unaltered is about -0.25 
eV. The remaining -0.28 eV are associated with config- 
urations with a double occupation of the N 2p "hole" 
orbital. Worthwhile to mention, the inter-pair corre- 
lation contribution involving the N 2p hole orbital and 
the remaining, doubly occupied N 2s, 2p orbitals arises 
from external two-particle excitations, i.e., excited dou- 
bly external configurations. This contribution amounts 
to about -0.84 eV. 

The findings above indicate that the total correlation 
energy of -4.00 eV contains a moderate contribution of 
about -0.28 eV originating from configurations in which 
the occupation of the hole orbital is no longer maintained. 
In order to obtain the loss of ground-state correlation, the 
difference of -3.72 eV should be compared with the cor- 
relation energy for the A'^-particle ground state. The cor- 
relation energy for the A'^-particle ground state is found 
to be about -4.77 eV. Hence, the correction to the di- 
agonal matrix elements, A7f^^(0), due to the loss of 
ground-state correlation shifts the N 2p bands downward 
by about 1.05 eV; see Table |l| 

The loss of ground-state correlation for the N 2s hole 
state is more difficult to calculate because maintaining 
the hole at the N 2s orbital in all configurations of the 
SDCI wave function is technically not possible. We cal- 
culated the associated correction by using only the corre- 
lation energy contribution of the excited doubly external 
configurations in the wave function I^'r"™^.') of the N 
2s hole state. The energy contribution arising from ex- 
cited singly external configurations is predominantly due 
to semi-external two-particle excitations resulting into a 
double occupation of the N 2s orbital. Hence, this con- 
tribution is beyond the loss of ground-state correlation. 
With the current choice of the cluster, no internal config- 
urations contribute to the net correlation energy of the 
N 2s hole state. Our result is that the overall correlation 
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correction shifts downward the center of gravity of the N 
2s bands by about 1.08 eV. 

Similar analyses were also carried out for the 
conduction-band (A^+l) states. The cluster exploited 
in these studies is [B13N28], with an active region [BN4]. 
The basis sets employed are 6-31 IG* for the B atoms and 
7-311G* for N. The SDCI wave functions for the (A^-l-l) 
electron states are constructed by correlating explicitly 
the 2s and 2p orbitals of the four N atoms around the 
central B site. The reference wave functions in the corre- 
lation treatment are the relaxed {N + 1) wave functions 

Given the configuration selection in MOLPRO, in a 
quasiparticle picture, the correlation between the doubly 
occupied N 2s and 2p orbitals and the singly occupied B 
2s-like or 2p-like orbital can be accounted for by SDCI 
wave functions where only the doubly occupied N 2s and 
2p orbitals are correlated. The singly occupied B 2s-like 
or 2p-like orbital is placed in the active orbital space, but 
no excitations out of this orbital are included in the cor- 
relation treatment. Such a SDCI wave function incorpo- 
rates correlation contributions of all relevant internal, ex- 
cited singly external, and doubly external configurations. 
The selected excitation space contains thus internal exci- 
tations and semi-external two-particle excitations which 
promote an electron from a doubly occupied N 2s or 2p 
orbital to the B 2s or 2p orbital and a second electron 
from the doubly occupied to the virtual orbital space. 
External one-particle and two-particle excitations which 
involve electrons or electron pairs from the doubly occu- 
pied N 2s and 2p orbitals are also included in the excita- 
tion domain. In analogy with the considerations for the 
(A''— 1) states, we compared the total correlation energy 
for the (A^-t-1) state, corresponding to SDCI wave func- 
tions of the type described above, with the correlation 
energy for the TV-particle ground state. The SDCI wave 
function for the iV-particle configuration is constructed 
by correlating the doubly occupied N 2s and 2p orbitals of 
the four N atoms around the central B site. The vacant B 
2s-like and 2p-like orbitals are all part of the variational 
orbital space in this case. 

The corrections to the diagonal matrix elements, 
AHl^^'^{0), due to the loss of ground-state correlation 
are similar for the B 2s and 2p electron-addition states, 
about 0.5 eV. The results are listed in Table HIl 

Although beyond the scope of this study, we mention 
that SDCI wave functions containing configurations in 
which the added electron is not maintained in the same 
B 2s or 2p orbital as in the reference wave function are 
related to satellite structures in the inverse photoioniza- 
tion spectra. These configurations describe correlation 
effects beyond the loss of ground-state correlation. Such 
SDCI wave functions can be constructed by including 
in the correlation treatment not only the doubly occu- 
pied N 2s and 2p orbitals, but also the singly occupied 
B 2s-like or 2p-like orbital. The inclusion of the singly 
occupied B 2s or 2p orbital in the correlation treatment 
leads to an additional contribution of about -1.0 eV. The 



latter effect is predominantly due to excited doubly ex- 
ternal configurations involving the extra electron. Such 
differential correlation effects lead thus to a stabilization 
of the (iV-|-l) electron configuration with respect to the 
iV-particle ground state by -0.5 eV. 

At this point, we have computed all relevant 
correlation-induced corrections to the diagonal matrix el- 
ements for the valence- and conduction-band states. The 
overall effect is a reduction of the HF band gap of c-BN 
from 13.62 eV to 6.91 eV. This result c ompa res well with 
experimental estimates of 6.0±0.5 eV^ZMl ^nd 6.4±0.5 
eV.ESI More sophisticated basis sets are not expected to 
affect significantly the calculated values, since the present 
choice of basis sets proves to be quite reasonable. For 
comparison, LDA calculations carried out with the triple- 
zeta GTO basis set used in the HF calculations yield a 
too small band gap of 4.34 eV, in agreement with other 
LDA studies; see, e.g., Refs. I70I71I An exception makes 
the LDA+GW study from Ref. |74| which provides a 
theoretical value of 6.3 eV for the indirect band gap of 
c-BN. 



B. Off-diagonal matrix elements of the local 
Hamiltonian 



We discuss next the correlation-induced corrections to 
the band widths. The computation of these corrections 
requires the off-diagonal matrix elements of the effective 
Hamiltonian in Eq. (10 1. 



having the hole at distinct sites Rj, see Eqs. (|5l^ 



We consider first the Hamiltonian matrix elements be- 
tween two frozen-orbital ROHF wave functions |<i>D~i, ;) 

^ The 

off-diagonal Hamiltonian matrix elements ^^rJ be- 
tween such mutually orthogonal wave functions l^^J^a') 
and \^o^i]ji) constitute the so-called hopping terms in an 
orthogonal tight-binding approach. The same considera- 
tions hold for the wave functions \^^^^^) and Hamilto- 
nian matrix elements H^^^^, . 

In a second step, we incorporate relaxation and po- 
larization effects in the nearby surroundings of the hole 
or extra electron. This is achieved by separate SCF op- 
timizations for each of the (A^— 1) and (A^-l-l) electron 
states. Hence, we construct the correlated wave func- 



tions |<i> 



N-l 



Af-l 
OmV 



and |<i> 



N+l 



l<i> 



On'al 



The sep- 



arate SCF optimizations lead to sets of non-orthogonal 
orbitals. A scheme for the computation of the Hamil- 
tonian and overlap matrix elements between such non- 
orthogonal wave functions was recently implemented in 
MOLPRO by Mitrushchenkov and Werner.'^' It is based 
on non-unitary transformations of the initial set of 
non-orthogonal orbitals to biorthogonal sets and follows 
closely the approach proposed by Malmqvist.^^ Alterna- 
tive approaches were proposed in other groups, see, e.g., 
Ref.. 98 

There are two different routes into constructing the 
correlated energy bands ek'^o-' and eki/o-. The first pos- 
sibility is to diagonalize a k-dependent matrix like that 
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in Eq. (10 1, which contains the Hamiltonian and overlap 



matrix elements between the non-orthogonal, correlated 
wave functions ) or )• A second route rests 



on deriving from the initial inter-site Hamiltonian and 
overlap matrix elements a set of effective hopping inte- 
grals for an orthogonal tight-binding-like formulation. In 
the latter approach, the effective hopping terms can be 
directly compared with the HF hoppings, which offer a 
more intuitive picture onto the effect of correlations on 
the electronic band structure. 

For electron-removal states which have different bind- 
ing energies and are non-orthogonal, i.e., iJo.m'm' 7^ 
H-Rj^mm and S'rj, mm' 0, the effective hopping matrix 
elements are defined as 



^mm'(T^j) 



1 



1 - S'2 

^ "^Rj, mm' 



R,j ,m7n' 



c2 

*~^R,j , mm' 



Ho .m'm' 



^ ) A ip2 {Huj^mm ~ -ffo.m'm')^ 1 
^ [ ^ -^Rj, mm' J 

and analogously for i„„'(Ri). Here, AE is the energy 
separation between the two eigenstates of the 2 x 2 CI 
secular problem, \H—ES\ — 0. This type of secular prob- 
lem in terms of non-orthogonal sets of orbitals will be re- 
ferred to as non-orthogonal CI (NOCI); see, e.g.,.'^If the 



two wave functions arc orthogonal, ($q, 
the effective hopping term reduces to 



N-l 

im'o"' 



N-l 
Rjmcr' 



!=0, 



'(Rj) 




AE — {Hflj^mm — Ho^m'rn'Y 



To compute the nn and nnn effective hoppings associ- 
ated with the N 2s, 2p valence bands, we designed two 
different clusters, [N20B40] and [N26B50]. The active re- 
gions of these clusters are the [N2B7] and [N3B10] frag- 
ments, respectively. The [N2B7] fragment consists of two 
N sites, denoted as 1 and 2 in Fig. [8] plus their nn B 
atoms. Likewise, the [N3B10] active region contains the 
N sites 1 and 2 plus the N atom denoted as 3 in Fig. |8] 
The nn B atoms of the three N sites are also included in 
the active subunit. The buffer regions of the two clusters 
consist of all N and B atoms in the first two coordination 
shells of the active boron sites. The [N26B50] cluster is 
sketched in Fig. [8] 

The basis sets exploited in the calculations for the 
[N20B40] cluster are 7-31 IG* for the N atoms and 6- 
311G* for B. In the case of the larger [N26B50] cluster, we 
used the 6-41G basis set for the boron atoms since oth- 
erwise the computation becomes unfeasible. In analogy 
with the analysis for the diagonal matrix elements, basis 
set effects for the effective hopping terms were checked 
by additional calculations for the nn hoppings using the 




FIG. 8: Illustration of the [N26B50] embedded cluster used for 
the computation of the nnn effective hopping terms associated 
with the valence-band N 2s, 2p hole states. The active region 
of the cluster consists of the [N3B10] fragment. The three N 
sites in the active region are depicted as large black spheres 
and labeled as 1, 2, 3. N atoms in the buffer zone are shown 
as small black spheres. The ten B atoms in the active region 
are drawn as large white spheres, whereas the buffer B atoms 
are shown as smaller white spheres. 



[N20B40] cluster. In this extra set of calculations, the 
smaller B 6-4 IG basis set was used. We found negligibly 
small basis set effects for the nn hopping terms, of the 
order of 1%. 

An overview of the valence-band nn and nnn effective 
hoppings is provided in Table |III| We list in the first 
column the Wannier functions where the hole resides. In 
the second column, we include results of frozen-orbital CI 
calculations in terms of the initial HF WO's (FOCI). In a 
next step, we allow for full relaxation of the doubly occu- 
pied 2s and 2p orbitals of the N atoms within the active 
regions (RO-NOCI). The orbital optimizations were car- 
ried out separately for each electron-removal wave func- 
tion I^Sfri,,) and yielding the wave functions 



I^Rj mo-'/ 



Rjmcr'/ 



and 1$^-/., 



Oni' a' 

respectively. The corresponding 
effective hoppings are listed in the third column. Finally, 
results where in addition to the relaxation of the N 2s, 2p 
doubly occupied orbitals, we also relax the hole orbital 



are given in the fourth column of Table [Hi] (ROH-NOCI). 

The results indicate small changes of about 0.02 eV for 
the nn effective hopping terms when relaxation effects in 
the nearby surroundings are accounted for. The general 
trend is a slight increase of the nn hoppings tmm' , which 
persists further upon the relaxation of the hole orbital. 
In the latter case, the changes are also of the order of 
0.01 eV. 
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TABLE III: Nearest-neighbor (nn), Rj = (0,1, -l)f, and 
next-nearest-neighbor (nnn), Rj = (0,1, 0)a, effective hop- 
ping matrix elements for the valence- band states (in eV). 
Frozen-orbital CI results (FOCI ) are given in the second col- 
umn. NOCI matrix elements in terms of relaxed orbitals are 
listed in the third column (RO-NOCI). NOCI results for which 
the orbital where the hole resides is also relaxed are listed in 
the fourth column. Matrix elements smaller than 0.01 eV are 
not included in the table. 



Active WO's 




nn tmm'(Rj) 






FOCI 


RO-NOCI 


ROH-NOCI 


2s-2s 


0.520 


0.541 


0.544 




0.184 


0.183 


0.187 


2Py(z)-'2Py(^) 


0.724 


0.739 


0.747 


2py(^y2p^(y) 


0.928 


0.954 


0.967 


'^Px{y)-2Py(x) 


0.053 (0.063) 


0.060 (0.070) 


0.062 (0.072) 


2Pa:(z)-2Pz(j;) 


0.053 (0.063) 


0.060 (0.070) 


0.062 (0.072) 


Active WO's 




nnn t„™/(Rj) 






FOCI 


RO-NOCI 


ROH-NOCI 


2s-2s 


0.011 


0.012 


0.012 


2px-2px 


0.054 


0.050 


0.049 


2py-2py 


0.037 


0.030 


0.028 


'2Py(z)-2pz(y) 


0.053 


0.050 


0.072 



Owing to the localized nature of the valence-band hole 
states, the effective hopping matrix elements decay fast 
with the distance Rj. The nnn hopping terms are found 
to be an order of magnitude smaller than the nn matrix 
elements. The short-range relaxation and polarization 
corrections for the nnn terms are also an order of mag- 
nitude smaller as compared to those for the nn elements, 
see Table Iml 

Some trial calculations were also performed to deter- 
mine the role of correlation effects beyond relaxation and 
polarization on the hopping matrix elements. As for the 
case of diagonal matrix elements, we employed CI wave 
functions with single and double excitations from the re- 
laxed 2s and 2p orbitals of the active N sites. We found 
that the SDCI treatment leads only to a slight reduction 
of the effective valence-band hoppings, with variations 
of up to 0.02 eV. Since these additional corrections are 
very small and the details regarding the construction of 
CI expansions where correlation effects beyond the loss 
of ground-state correlation are minimized require very 
technical explanations, we do not discuss these effects 
any further. 

Summarizing our results for the effect of correlations 
on the nn and nnn intersite interactions, we find that the 
widths of the upper-lying valence bands are only little af- 
fected. The fact that the hopping matrix elements from 
the calculations in terms of relaxed orbitals are slightly 
enlarged as compared to the FOCI values is due to the 
larger intersite overlap between the relaxed orbitals of 



TABLE IV: Nearest-neighbor (nn), Ri = (0,-1, l)f, and 
next-nearest-neighbor (nnn), Ri = (0,-1, 0)a, effective hop- 
ping matrix elements for the conduction-band states (in eV). 
FOCI results are listed in the second column. The third col- 
umn summarizes the NOCI data. Matrix elements smaller 
than 0.01 eV are not included in the table. 



/active VVVJ S 




nn tn„' (ri.1 ) 




FOCI 


RO-NOCT 




n Dzin 


n ofii 

U.UU 1 


^Px -^Px 


1 '^7^ 

1 .O 1 -L 


1 fi9fi 


9t) / \-9ti / \ 


0.383 


0.377 


2'D^( ^,\-2v^,i ^\ 

rxi^y) ^yy{^x) 


0.269 


0.255 


'2Px(z)-'2Pz(x) 


0.269 


0.255 


'2Py(z)-'2Pz(y) 


0.168 


0.190 


Active WO's 




nnn t„„/(Ri) 




FOCI 


RO-NOCI 


2px-2px 


0.177 


0.164 


2pz-2pz 


0.131 


0.116 


2py-2Py 


0.070 


0.063 


2Py(x)-'2Px(y) 


0.031 


0.030 


2Px(z)-'2Pz(x) 


0.062 


0.060 



hole resides are polarized toward the positive charge at 
this site. It is this orbital polarization effect that leads 
to a larger intersite orbital overlap and enhanced hop- 
pings. Similar polarization effects were also discussed for 
the case of O hole states in MgO.I^ We have not con- 
sidered explicitly the effect of long-range polarization on 
the off-diagonal matrix elements. This effect, however, is 
expected to be negligible. 

A similar analysis was carried out for the lowest 
electron- addition states and conduction-band hoppings. 
The clusters employed to investigate the nn and nnn B- 
B interactions are designed in perfect analogy with the 
[N20B40] and [N26B50] clusters described above. The ac- 
tive regions of these [B20N40] and [B26N50] clusters are 
the kernels [B2N7] and [B3N10], respectively. Results of 
both FOCI and NOCI calculations are summarized in Ta- 



the two wave functions I'&RjJjg./) 



and |<i>, 



N-l 



The N 



2p orbitals within the nn region of the N site where the 



ble IV In the NOCI calculations, the 2s and 2p orbitals 
of the nearest N atoms are allowed to relax in response 
to the extra electron in the conduction-band B 2s or 2p 
orbital. The additional effects resulting from the relax- 
ation of the orbital where the extra electron resides are 
difficult to estimate for the nn hopping terms associated 
with the (A^-l- 1) states because the orbital could not be 
restrained from delocalizing over the two B sites in the 
{N + 1) SCF calculation. Such effects could be assessed, 
however, for the nnn hoppings t„„' and proved to be of 
the same order of magnitude as that observed for the 
valence-band (N — l) states, i.e., of the order of 0.01 eV. 

The larger extent of the conduction-band WO's mani- 
fests itself in the larger values of the nn and nnn effective 
hopping terms and their somewhat slower decay with the 
distance Rj. The orbital relaxation and short-range po- 
larization effects on these matrix elements, though, re- 
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main relatively small, with corrections of about 0.02 eV. 
The largest change is observed for the 2px-'2,px intersite 
interaction, about 0.06 eV. Our results indicate modi- 
fications of the hopping terms in two directions: some 
matrix elements are slightly enlarged by taking into ac- 
count the short-range relaxation and polarization, others 
are reduced. Similar findings were reported for the effec- 
tive hoppings associated with the lower-lying Mg 3s and 
3p conduction bands of yigO^^ 

For the computation of the nnn hopping terms, the 
basis set employed for the B atoms is the smaller 6-41 G 
set. Basis set effects were assessed by extra calculations 
for the nn hoppings, using either 6-41G or 6-31 IG* basis 
sets for the boron species. The cluster employed in these 
calculations is the [B20N40] cluster. With the 6-41G ba- 
sis set, we find changes of 0.01 to 0.06 eV for the nn B 
2p-2p hoppings and changes of about 0.05 eV for the nn 
B 2s-2s hoppings, both at the frozen-orbital HF and cor- 
related levels. Basis set effects of similar magnitude, not 
larger than few tens of meV, are expected for the nnn 
conduction-band hoppings. 

To conclude, the overall corrections for the intersite 
interactions, related to the formation of the short-range 
polarization cloud around the added electron, are rather 
small, in the range of few percents. A similar situation 
has been also observed for the conduction-band states 
of the ionic insulator MgO.f^^The widths of the HF va- 
lence and conduction bands change thus little by incor- 
porating short-range relaxation and polarization effects. 
In contrast with the relatively weakly correlated BN 
compound, the effective hopping matrix elements associ- 
ated with electron-removal and electron-addition states 
in strongly correlated transition-metal compounds like 
the perovskite manganites are reduced considerably, by 
a factor of 4, due to the presence of strong magnetic in- 
teractions in these systems A strong reduction of the 
effective hopping terms was also found for the case of 
cuprate compounds .^221 

At last, we evaluate the effect of the correlation- 
induced corrections to the widths of the bands on the 
size of the band gap. For this purpose, we consider first 
the energy of the N 2p band states at the F^g point. 

For the Px band, for example, eo^,i=£'i:'?S°'' + Sii^x*^-* + 
4ti°^^'' -I- 0{nnn) + ... where £'1*^°"' is a diagonal matrix 
element. For symmetry reasons, t^x]c^^ = t^,l^^ = tf,l^\ 
Taking into account corrections of about 0.02 eV for ii^i^'' 
and 0.003 eV for t^°l^\ we find that the N 2p bands are 
shifted upwards by about 0.17 eV at the F point. The 
band gap is reduced then by the same amount. For the B 
2s, 2p conduction-band complex, the correlation-induced 
corrections to the B 2p-2p intersite matrix elements are 
most relevant at the Xic symmetry point. The four by 
four secular problem factors out on account of symmetry 
at the X point, such that one finds the degenerate B 2py 
and 2p2 states and the lower 2s and 2px bands. The en- 
ergy of the 2px band at the symmetry point Xic (k^ — 

0, ky = TT, fc^ = tt) is given by ekAt,2^=£^i°°°'' + 4ii°"^ - 



St^^^ - 2ii°i°) -2^^ + 2tp^ + ....IMlFrom symmetry 
considerations, ti]l°'' = 4°"'' = 4°i"^\ ii™"^^ = and 
t'"x!x'^ = t^,y^^ ■ Corrections of about 0.05 eV for t^^x-^"^ and 
0.01 eV for and for the nnn hopping terms lead to 

a shift to lower energy of the B 2p band at the Xic point 
and a further reduction of the gap by about 0.24 eV. Our 
final result for the fundamental gap of BN is thus 6.50 
eV, which compares well with experimental estimates of 
6.0±0.5 eV-2^and 6.4±0.5 eV.I^ 



V. CONCLUSIONS 



We investigated in the present work the correlation- 
induced corrections to the Hartree-Fock band structure 
of zinc-blende BN. A local Hamiltonian approach and 
quasiparticle approximation were employed for this pur- 
pose. An important advantage of our wave-function- 
based method is that it allows for a rigorous, ab initio 
determination of the relevant correlation contributions to 
the gap and band widths of a crystalline material without 
compromising its infinite nature. The correlation contri- 
butions associated with the short-range relaxation and 
polarization around an extra hole or electron were ex- 
plicitly computed. Moreover, correlation effects beyond 
relaxation and polarization, such as the loss of ground- 
state correlation, were also investigated. The long-range 
polarization of the crystal was treated within the dielec- 
tric continuum approximation. 

We found that the on-site and nearest-neighbor relax- 
ation and polarization bring contributions of few eV to 
the correlation-induced corrections to the gap of BN. The 
long-range polarization caused by the extra hole or elec- 
tron leads to an additional reduction of about 2.0 eV of 
the HF band gap. The net effect is a reduction of the HF 
gap from 13.62 eV to a value of 5.36 eV. If, in addition, 
we take into account correlation effects beyond relaxation 
and polarization, such as the loss of ground-state correla- 
tion, the band gap increases to 6.91 eV. At last, an addi- 
tional small reduction of 0.41 eV to a value of 6.50 eV is 
caused by correlation-induced variations of the widths of 
the bands. The computed band gap compares well with 
that deduced from soft x-ray experiments at the B and N 
iC-edge^^^EII and optical absorption measurements.!^ As 
expected, the LDA-DFT calculations yield a too small 
band gap of 4.34 eV. 

While the HF gap is substantially modified by cor- 
relations, the widths of the bands are less sensitive to 
such effects. The inclusion of local correlations leads to 
a slight broadening of the N 2p bands. This broadening 
is attributed to the polarization of the nearest-neighbor 
N 2p orbitals toward the site where the hole resides. 
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APPENDIX: 6-311G* AND 6-41G BASIS SETS 
FOR THE B ATOMS. 

Tables IVl and I VII list the two basis sets for the B atoms 
designed for the computation of the correlation-induced 
corrections to the HF gap of c-BN. 



TABLE VI: Exponents (in a.u.) and coefficients of the B 
6-41G GTO basis set. 



Shell 


Type 


Exponent 


Coefficient 
s p 


2 


sp 


13.9313977 


0.0853578 


0.0418012 






3.1396444 


0.7038039 


0.2468341 






0.9037473 


0.7225717 


0.7649465 






0.4893080 


-3.1330026 


0.7122470 


3 


sp 


0.1691601 


1.0 


1.0 



TABLE V: Exponents (in a.u.) and coefficients of the 6-311G* 
GTO basis set used for the B atoms. 



SheU 


Type 


Exponent 


Coefficient 
s p 


2 


sp 


13.5055949 


0.1099474 


0.0420302 






3.0900933 


0.8728330 


0.2359840 






0.9128560 


0.5149721 


0.7318254 


3 


sp 


0.4927046 


1.0 


1.0 


4 


sp 


0.213 


1.0 


1.0 
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The wave functions l^'^ ^ ,) and l^'iY^^ ) are correlated 
only through one-particle excitations. 
For each of these p-like lobes, two of the nitrogen atoms 
of a BN4 tetrahedron are closer as compared to the other 
two N sites, see Fig. [S] 

Removing the polarization d function has negligible effects 
(<<1%) on the on-site, nn, and rmn orbital relaxation 
contributions. 

\AHI^^'^{0)\ increases by only 2 % when enlarging the 
[N55B28] cluster to [N63B44]. 



